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Abstract 

We present new family of exact analytic solutions for three anyons in a 
harmonic potential (or in free space) in terms of generalized harmonics on 
5^, which supplement the known solutions. The new solutions satisfy the 
hard-core condition when a = |, 1 (a being the statistical parameter) but 
otherwise, have finite non-vanishing two-particle colliding probability den- 
sity, which is consistent with self-adjointness of the Hamiltonian. These so- 
lutions, however, do not have one-to-one mapping property between bosonic 
and fermionic spectra. 
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Smooth interpolation between bosonic and (spinless) fermionic spectra has been a very 
attractive idea and appeared the concept of anyon [|1| based on the homotopy group in two 
space dimensions. The anyons are described effectively by Schrodinger equation of bosonic or 
fermionic particles with statistical gauge field of Aharonov-Bohm type. One can also derive 
this equation of motion from the second quantized abelian Chern-Simons gauge theory. 
It is believed that fractional quantum Hall effect is an example of physical phenomenon of 
anyons P|. 

On the other hand, the detailed analysis of many anyons has not been successful. Even 
for simple system such as harmonic oscillator (HO), analytically known exact solutions 0-0] 
fall short of the numerical result in |^ which demonstarates smooth interpolation. Only 
perturbation approach could reproduce the numerical result. The main difficulty of the 
problem lies in the non-trivial exchange property of anyons and hard-core condition on its 
wavefunction. 

In this Letter, we present new family of analytically exact solutions of three particles 
in HO or in free space. These solutions supplement the presently known analytically exact 
solutions but unfortunately do not satisfy the hard-core condition except at a = 1/3, 1 (a is 
the statistical parameter chosen as < a < 1, a=0 at for boson and 1 for spinless fermion). 
None the less, considering the situation of absence of analytic solutions for missing states, 
we think these solutions will provide valuable information about anyons. 

Our analysis is done by adopting ordinary Hamiltonian with anyonic particle-exchange 
property imposed on wavefunction. We use coordinates system following Coordinates 
of three particles are represented as complex numbers, Za = Xa + iya where a = 1,2,3. 
We separate out center of mass (CM is irrelevant to the anyonic property) by using CM 
coordinate, Z = and relative ones u = (^i+^W^s) ^^^^ ^ ^ (ziW^+vz,) ^-^^^ 

i2TV 

rj = e 3 . 

Relative motion (RM) of free particles is described in terms of Hamiltonian (scaled 
as dimensionless) , Hfree = ~ considering systems with non-singular 

potential energy of the form, V = V{r) where r = y/uu + vv > 0. This includes the free 
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case, V — and HO case, V — r'^. This Hamiltonian is invariant under exchange of any two 
particles. To see this in our notation, let us denote the second and third particle exchange 
operation, (1, 2, 3) — > (1, 3, 2) as E, and cyclic operation, (1, 2, 3) — > (2, 3, 1) as P. Any 
two-particle exchange is represented in combination of £"s and P's. The definition of u and 
V shows that P : {u,v) {rfu,r]v) and E : {u,v) — > {v,u). Therefore, RM is trivially 
invariant under P and E. 

For separation of variables, we introduce spherical coordinates, u — rsin(^/2) e^^^+a) 
and V — rcos(^/2) e*^^""^\ ^, 9 and (f) are angle variables whose fundamental domain is 
given as < r, < ^ < |, — | <0<|, 0<6'<27r since under P and E, we have 
P : (r, ^, 0, ^ (r, ^, + ^, ^ + tt) and E : (r, ^, 0, 6) ^ (r, tt - ^, -0, respectively. The 
angular domain is not but ^Xz-^ ■ (Euler angles on 5"^ or on 5" C/( 2) are defined typically as 
(C) X = ~4>j ip = 29 — 2tt) whose ranges are given as0<(^<7r, 0<x< 27r, — 27r <ip< 2tt). 

Anyon wavefunction is defined to have a phase e*'^'^ when any of two particles are inter- 
changed. This requires the wavefunction have the phase under P and E, 

E: ^'(r,7r-e,-0,^) = e-'^*(r,e,0,^) 
27r 

P : *(r, e, + — , ^ + tt) = (^^^-^{r, <P, 0) 

R : *(r, e, 0, ^ + 27r) = e^6'^'^*(r, 0, 9) (1) 

The operation i?, 27r rotation of ^ corresponds to a composite operation EPEP. 
The Hamiltonian is written in (r, ^,4>,9) coordinates as 

"—itr.^^'-'^^^ir.M^vir^) (2) 

where M is a Laplacian on S^, 

^d^^'^^^d^^ sm^{C/2y2id9^1W ^cos2(e/2)^2i9^ 7^0^ ' ^ ^ 
M commutes with the Hamiltonian H. The relative angular momentum, L, 



is another useful mutually commuting operator. A simultaneous eigenstate of H , M and 
L is given in factorized form, \E'£;_^_i(r, 0, 6') = RE,^iij') S^.K'C; 0? is a harmonic on 

S*^, whose two-particle analogue is e^*°^^^ where < ^12 < tt is the angle of the relative 
coordinate. E, yu(/i + 2) and I are eigenvalues of if, M and L respectively. 

Since one can easily solve the radial solution once the harmonics are found (for example, 
it is written as a Bessel function of the first kind for free particle and as a Laguerre for HO), 
we concentrate on finding the angular part. Boundary conditions for S are obtained ||7|JT0(| 
from Eq. (||), consistently with self-adjointness of L and M, 

S(e,</),^ = 27r) = e^''"^ S(e,<^,e = 0) 

= ^, = vr) = ^ = ^ = q) 

|lnS(e,0^^,.^.)^^lnH(e,0^^,.^O) 

I lnS(^ = |,0,0) = -| lnS(e = |,-0,^^), (5) 

where last two identities hold for < </> < vr/S. In addition, any current across the boundary 
needs to be finite, which requests 

finite. We note that even though hard-core condition (vanishing at ^ = t{/2 and = 0) is 
consistent with this boundary conditions, it is not the unique choice for a 7^ odd integer as 
we will see promptly. 

The harmonics are not single- or double-valued but multi-valued since boundary condi- 
tions contain e*""^ instead of ordinary ±1. To find the harmonics, we employ ladder operators 



TT| satisfying the su{2) Lie-algebra, [K^,K^ = ±K± and [Kj^,K_] = 2K^ where 



K c^(>\i^ ^ ^ cot^g 

+ ^ di sine 90 2 dQ^ 

na r d Id cot f (9 , 

(74 sm 4 (70 2 o9 
1 (9 1 ^ 



where = K+. (Integration measure is sin 2^d^d(f)d6). K± and are invariant under 
P and R but get extra - sign under E. Noting that M = 4{ K+k^+K-K+ ^ x^}, /i/2 is 
Casimir number for harmonics and is semi-positive definite. 

Well-known analytic solutions called type / (//) are reproduced in our approach 
if we apply K"^ (-^+) on states which are annihilated by (-^-)- Of course, they satisfy 
the hard-core condition for < a < 1 and fi = integer +3a > for type I {fi = integer 
—3a > 3 for type II). 

Newly found solutions cannot be constructed this way. They do not have states which 
are annihilated by or K^. Suggestion of this new type has been made in |T^. Similar 
solutions but approximate ones were considered in These have /i = integer ±a > 2 
(-|- (— ) for type /// {IV)). (Note that the convention for the class of the solution in |^ is 
different from ours: //, /// and IV corresponding to our ///, IV and II respectively) . 

A solution with (/x = 3 — a, / = — 3-|-3a) (type IV), resulting in the interpolation between 
the fermionic ground state and a bosonic excited state of HO, is given as (using notations 
of z = v/r and w = u/r) E^^a-s+sa = {z^ - u^^)"-P3-q -3+3a with 

_|_ 1^)3-20 

Ps-a _3+3„ = ■ — — + 2 more cychc permutations , (8) 

[Z^ + zw + w^)" 

where 2 more added terms are obtained from the first term if z {w) is replaced by rjz 
ijfw), and by rfz {rjw) respectively. Note that for \w\ < \z\, rfw = w but rj^z ^ z 
because of multi-valuedness of the factors in 'B.-j^-a-z+za- The solution is equally put as 
S3_a _3+3a = {z — w)'^{z + + {z — i]w)°'{z + r/tZ;)^~^" + {z — rf'w)'^{z + ffw)'^''^'^ and 

more symbohcally as S3_q,^_3+3q, = S{{z — w)"^{z + Its explicit form is given as 

S3^_3 = 3(z^ -|- w^) in the bosonic limit, and ^2,0 = 3{zz — ww) in the fermionic limit. 

One can check that H3_q, „3+3„ is an eigenfunction of M using Kj^z = —iw, K^w = iz, 
K_z = —iw and K_w = iz. E^^a-s+sa satisfies the boundary conditions in Eqs. |^), yet 
does not vanish but is finite when two of three particles coincide {z = w = l), except for 3a = 
odd integer: 

E,.^,.s+3a{z = w = 1) = -(1 - r^)" (1 + e^^™) . (9) 



We have used identities such as hm^^oll + (1 ~ = ^'5, hm5_,o(l + (1 ~ ^)v'^) = ^^'^ 
and their complex conjugates. (Coordinates {zi, Z2, z^) is not suited to proving Eq. @ 
due to lack of faithful representation of (—1)"). Two particles in the presence of the third 
can collide each other, as seen in two-particle case in the context of self-adjoint extension 
||13|| . The chance for simultaneous collision of three particles is, however, strictly prohibited 
except the bosonic limit because radial part vanishes. 

Will the enumeration of these solutions (FIG. 1, whose details will appear elsewhere) 
lead to smooth interpolation? A solution is missing, interpolating between (/i=2, /=— 2) 
boson and (/i=3, /=1) fermion, but this only implies that fi^2+a. A real potential problem 
against the smooth interpolation comes from the possibility of a two-to-one correspondence. 
For example, two bosonic states with {fj.=5, l=—l) and (/i=3, /=— 1) are mapped into a 
fermionic state with (/i=4, 1=2): 

H5-«,-i+3„ = S{iz- wf+'^iz + , (10) 

S3+«,_i+3a = S{{z- wY^^iz + wY+'-} . (11) 

S5_Q _i+3a is an extra solution which has no counter-part to the numerical solution in Ref 
ip. It satisfies the boundary conditions in Eqs. (^ P) and is obtained by applying K"^ on 
H5_q„5+3q, = S {{z — w)" {z -\- wY~'^°'} . Since that hard-core condition is satisfied at a=l/3, 
the relaxed boundary condition cannot be the cause of this trouble. It might be possible 
that the extra solution is an isolated one which cannot be interpolated smoothly to bosonic 
and fermionic limit satisfying the hard-core condition. Even so, it does not eliminate the 
fact that there are extra states beyond smooth interpolation. 

For HO, energies of our solutions are linear in a, since the energy is given as /i + 2 + 2n 
{n is a non-negative integer and 2n is due to the radial part). However, as we impose the 
hard-core condition, the wavefunctions will adjust themselves and energies will have a slight 
non-linear dependence on a. This non-linear dependency should be manifest at the near side 
of the fermionic states of new solutions since all the solutions satisfy the hard-core condition 
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at a = 1/3. In fact, this property is seen in numerical 0, perturbative and approximate 
estimation p, even though /i at a = 1/3 does not agree completely. 

One can check that reflection symmetry in a by 2 — a for type / and //, and a 
by 1 — a for /// and IV ^ although the latter is not perfect, being subject to hard-core 
boundary condition. We also expect that the tendency of linear dependence of /i on a for 
anyons, ( ^^^~"^'' )q;, — 2)q;, ■ ■ -, — 2)q:, and _( ^(^~i) ^q, to an integer, will 

persist if we relax the hard-core boundary condition. This trend is shown in the numerical 
analysis for four anyons fTsH . 

In summary, we present four types of analytically exact normalizable wavefunctions for 
three anyons, two of which are newly found. The solutions are given in terms of generalized 
harmonics on S^, with boundary conditions consistent with anyonic exchange property, self- 
adjointness of the Hamiltonian and no-infinite current across the boundary. The hard-core 
condition turns out to be an additional constraint on the harmonics, which is obeyed by 
new family of solutions at a=l/3, 1. The analytic expressions for the whole solutions do 
not give a one-to-one correspondence mapping between the bosonic and fermionic spectra, 
whose way out is not seen at present. 

The author has a great benefit from discussions with Professors R. Jackiw, M. Rocek, J. 
Verbaarschot and I. Zahed, and Drs. G. Amelino-Camelia and K. Cho. 
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FIGURES 



(fli) o\ bosonic states (//J ) of fermionic states 




o^O a=l/3 a=l 



FIG. 1. Four types of interpolation between bosonic and fermionic states: Schematic diagram 
for Casimir number (/x) v.s. statistical parameter (a). Solid lines represent Type / and //, and 
dashed ones /// and IV. Wiggly lines are extras which prevent smooth interpolation. Exponents 
on brackets denote degeneracy. 
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